Accepted for publication in The Astrophysical Journal 

Preprint typeset using I^'T^]X style emulateapj v. 10/09/06 



O 

o 

< 

(N 
< 

o 

6 



> 

m 

m 

^' 

o 
o\ 
o 



EFFECTS OF PRIMORDIAL MASS SEGREGATION ON THE DYNAMICAL EVOLUTION OF STAR 

CLUSTERS 

Enrico Vesperini"'^, Stephen L.W. McMillan""^, Simon Portegies Zwart ^'^ 

^ Department of Physics, Drexel University, Philadelphia, PA 19104; vesperin@physics.drexel.edu, steve@physics.drexel.edu 
^ Astronomical Institute "Anton Pannekoek" and Section Computational Science, University of Amsterdam, Kruislaan 403, Amsterdam, 

The Netherlands. 

3 Sterrewacht Leiden, P.O. Box 9513, 2300 RA Leiden, The Netherlands; S.F.PortegiesZwart@uva.nl 
Accepted for publication in The Astrophysical Journal 

ABSTRACT 

In this paper we use TV-body simulations to study the effects of primordial mass segregation on the 
early and long-term evolution of star clusters. Our simulations show that in segregated clusters early 
mass loss due to stellar evolution triggers a stronger expansion than for unsegregated clusters. Tidally 
limited, strongly segregated clusters may dissolve rapidly as a consequence of this early expansion, 
while segregated clusters initially underfilling their Roche lobe can survive the early expansion and 
have a lifetime similar to that of unsegregated clusters. Long-lived initially segregated clusters tend 
to have looser structure and reach core collapse later in their evolution than initially unsegregated 
clusters. We have also compared the effects of dynamical evolution on the global stellar mass function 
(MF) of low-mass main sequence stars. In all cases the MF flattens as the cluster loses stars. The 
amount of MF flattening induced by a given amount of mass loss in a rapidly dissolving initially 
segregated cluster is less than for an unsegregated cluster. The evolution of the MF of a long-lived 
segregated cluster, on the other hand, is very similar to that of an initially unsegregated cluster. 
Subject headings: globular clusters: general, methods: n-body simulations, stellar dynamics 



1. INTRODUCTION 

Mass segregation, the tendency of more massive stars 
to preferentially populate the inner parts of a star clus- 
ter, is one of the consequences of two-body relaxation 
and of the evolution toward energy equipartition in stel- 
lar systems. The characteristic timescale for this process, 
Tms, is, for a population of massive stars with mass m^, 
Tms ~ {{m) / ruhjUeiax where Ueiax is the cluster relax- 
ation time and (m) the mean mass of the cluster stars. 

However, a number of young clusters with ages sig- 
nificantly smaller than the time needed to produce the 
observed mass segregation by standard two-body relax- 
ation show a significant degree of mass segregation (e.g. 
Hillenbrand 1997, Hillenbrand & Hartmann 1998, Fis- 
cher et al. 1998, de Grijs et al. 2002, Sirianni et al. 2002, 
Gouliermis et al. 2004, Stolte et al. 2006, Sabbi et al. 
2008). 

The results of a number of theoretical studies (e.g. 
Klessen 2001, Bonnell et al. 2001, Bonnell & Bate 2006 
but see also Krumholz et al. 2005, Krumholz & Bon- 
nell 2007) showing that massive star would preferentially 
form in the center of star-forming regions suggest that 
the observed segregation in young clusters would be pri- 
mordial and imprinted in a cluster by the star formation 
processes. Possible dynamical routes leading to early 
mass segregation in young star clusters have been studied 
by McMillan, Vesperini & Portegies Zwart (2007). 

In this paper we present the results of a survey of A^- 
body simulations aimed at exploring the implications of 
primordial mass segregation for the dynamical evolution 
of star clusters. We have explored the evolution of clus- 
ters located at different galactocentric distances and with 
different degrees of initial mass segregation and com- 
pared the evolution of segregated clusters with that of 
clusters with the same initial density profile but no ini- 



tial mass segregation. Our models focus on the effect of 
primordial mass segregation on a purely stellar system; 
primordial gas and the effects of its expulsion (see e.g. 
Baumgardt & Kroupa 2007) are not included here. 

In Sect. [2] we present a preliminary analytical calcula- 
tion aimed at exploring the differences in the response of 
clusters with different initial degrees of mass segregation 
to the early mass loss due to stellar evolution. This is, in 
most cases, the first process affecting the cluster struc- 
ture well before the effects of two-body relaxation become 
important. We continue, in SectlSl with the presentation 
of the results of our iV-body simulations, focussing our 
attention on clusters lifetime, structural evolution and 
the evolution of the stellar mass function. We discuss 
and summarize our results in SectUl 



2. EARLY EVOLUTION OF SEGREGATED CLUSTERS: 
ANALYTICAL ESTIMATES 

A number of studies (see e.g. Chernoff & Shapiro 
1987, ChernoflF & Weinberg 1990, Fukushige & Heggie 
1995, Takahashi & Portegies Zwart 2000) have shown 
that early mass loss due to stellar evolution can have 
a significant impact on the evolution of clusters: this 
early mass loss causes a cluster to expand and, for a low- 
concentration cluster, leads to the cluster's complete and 
rapid dissolution. For initially mass-segregated clusters, 
the mass lost due to the evolution of massive stars is re- 
moved preferentially from the cluster's inner regions, and 
the early expansion of the cluster is stronger and poten- 
tially more destructive than when the same amount of 
mass is lost in a non-segregated cluster. 

Fig. 1 shows the results of a simple semi-analytical cal- 
culation illustrating the augmented destructive effect of 
stellar mass loss in a mass-segregated cluster. In this fig- 
ure we plot the virial ratio, T/V, of an isolated cluster 
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with an initial Plummer density profile 
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Fig. 1. — Virial ratio of a cluster with a Plummer density profile 
after the rapid loss of a AA/ fraction of its initial mass, versus 
AM. The density profile of the mass lost is assumed to be that of 
a Plummer model, but with a scale radius ajv/L equal to or less than 
the scale radius of the whole cluster, a^-i^gter- The number beside 
to each curve indicates the corresponding value of aML / o-cluster 



In order to mimic the preferential mass loss from the 
inner regions of a mass-segregated cluster, we have sim- 
ply assumed that the density profile of the mass lost also 
follows a Plummer model, but with a scale radius {uml) 
smaller than the scale radius of the cluster (aduster)- The 
new potential {V) and kinetic (T) energies for the den- 
sity profile thus obtained are calculated. The curves in 
Fig. 1 show the resulting ratio T/V as a function of AM 
for different values of um l/ aduster- 

For ajv/ l = o-duster one recovers the well-known result 
that a system becomes unbound for AM > 0.5 (Hills 
1980; see Boily & Kroupa 2003 for a more detailed cal- 
culation showing that larger values of the mass loss are 
actually necessary to unbind a cluster). However, for 
a cluster preferentially losing mass from the inner re- 
gions, uml < OLciusteri wc sce that a significantly smaller 
amount of mass loss, innocuous for an unsegregated clus- 
ter, can lead to the dissolution of a mass-segregated clus- 
ter. 

This simple semi- analytical calculation underscores the 
potentially crucial implication of initial mass segregation 
for the evolution of star clusters; the results of the N- 
body simulations presented in the rest of this paper il- 
lustrate these consequences in more detail. 

3. DYNAMICAL EVOLUTION OF SEGREGATED 
CLUSTERS: Af-BODY SIMULATIONS 

3.1. Method and Initial Conditions 



The results presented in this section are based on 
a survey of direct iV-body simulations carried out 
using the starlab package (Portegies Zwart et al. 
2001; http://www.manybody.org) and accelerated by a 
GRAPE-6 special-purpose board (Makino et al. 2003). 
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Fig. 2. — Initial radial density profile (the radius is normalized 
to the core radius, r^) for the S (solid red line) and the U (dashed 
black line) clusters. 
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Fig. 3. — Initial radial profile of the stellar mean mass for the 
initially segregated systems investigated in our simulations: 5 runs 
(solid red line/filled dots), M run (dotted green line/open trian- 
gles), L run (blue dot-dashed line/asterisk dots). 



We have followed the evolution of clusters located at 
different galactocentric distances, Rg = 1,4,18,40 kpc, 
in a Galactic tidal field modeled as a Keplerian poten- 
tial determined by a point mass Mg equal to the total 
mass of the Galaxy inside the distance Rg. We have 
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studied both tidally truncated clusters and clusters un- 
derfilling their Roche lobe with a ratio R/Rt of the total 
radius, R to the tidal radius, Rt equal to 0.75 and 0.5. 
A Kroupa (Kroupa, Tout & Gilmore 1993) stellar ini- 
tial mass function with star masses between 0.1 and 100 
7710 was adopted and star masses were generated by the 
analytical function in Eq. 14 of Kroupa et al. (1993). 

For segregated systems, mass segregation was set up 
by first letting the cluster evolve without including the 
effects of stellar evolution until a given degree of segrega- 
tion was reached due to normal two-body relaxation. We 
emphasize that this is just the procedure we have used 
to generate a self-consistent initially segregated cluster 
and it is not meant to model any stage of cluster evo- 
lution. Alternative procedures to produce initially seg- 
regated clusters have recently been suggested (see e.g. 
Subr et al. 2008, Baumgardt et al. 2008). 

Fig. 2 and Fig. 3 show, respectively, the initial density 
profile and the initial radial profile of the star mean mass 
for the systems we have studied; in particular, the solid 
line in Fig. [3] shows the initial degree of mass segregation 
of the standard models on which we focus in this paper. 
The initial number of particles for the standard segre- 
gated {S) and unsegregated [U) runs is = 25K. Initial 
conditions for moderate and low initial mass segregation 
(the M-Rgl8 and L-Rgl8 runs, respectively) are obtained 
by extracting snapshots at earlier times of the same pre- 
liminary simulation mentioned above when the system 
has a lower degree of segregation and has lost a smaller 
number of particles; the initial number of particles for the 
M-Rgl8 and the L-Rgl8 runs are, respectively, A^ — 33K 
and A = 38K. 

In order to explore the dependence on A^ of the evolu- 
tion of clusters most affected by initial mass segregation, 
simulations with the same degree of initial segregation 
as S-Rgl8, M-Rgl8 and L-Rgl8 have been repeated with 
a larger number of particles (A^ = 70K, N = 103K and 
A^ = 120K respectively). 



Id. 


Rg(kpc) 


R/Rt 


Initial mass segregation 


S-Rgl 


1 


1 


strong 


S-Rg4 


4 


1 


strong 


S-Rgl8 


18 


1 


strong 


S-Rgl8R075 


18 


0.75 


strong 


S-Rgl8R05 


18 


0.5 


strong 


S-Rg40 


40 


1 


strong 


U-Rgl 


1 


1 


no segregation 


U-Rg4 


4 


1 


no segregation 


U-Rgl8 


18 


1 


no segregation 


U-Rgl8R075 


18 


0.75 


no segregation 


U-Rgl8R05 


18 


0.5 


no segregation 


U-Rg40 


40 


1 


no segregation 


M-Rgl8 


18 


1 


moderate 


L-Rgl8 


18 


1 


low 



TABLE 1 

Initial conditions of all simulations. The columns list (1) 
Id (Standard Segregated Run-5, Standard Unsegregated 
Run- {/,MoDERATE Segregation Run-M, Low Segregation 

RUN-I/),(2) GALACTOCENTRIC distance, Rg IN KPC, (3) RATIO OF 
the INITIAL CLUSTER RADIUS TO THE TIDAL RADIUS R/Rt, (4) 
DEGREE OF INITIAL MASS SEGREGATION 

We summarize in Table 1 all the initial conditions con- 
sidered, along with the identification used throughout 
the paper to refer to each simulation. 



3.2. Results: cluster evolution and lifetime 

The analytical calculation presented in Sect 12] showed 
that early impulsive mass loss associated with stellar evo- 
lution can have a stronger impact on initially segregated 
cluster than unsegregated clusters. 

For tidally truncated clusters with fixed initial mass, 
evolving in a host galaxy with constant circular veloc- 
ity, the mean cluster density pud oc M/Rf oc Mg/Rg 
decreases with the galactocentric distance Rg as pud oc 
l/i?g. The cluster dynamical time, defined using the 
mean density within the tidal radius, t^d oc 1/ ^pud thus 
increases linearly with Rg (note that this dynamical time 
is proportional to the circular orbital period of the cluster 
around the Galaxy at distance i?g). This scaling implies 
that, for tidally truncated clusters, the amount of im- 
pulsive mass loss due to stellar evolution increases with 

Rg. 

Fig. [5] shows the dependence of the cluster dissolution 
time, Tdissi (defined as the time when 1 percent of the 
initial mass is left in the cluster) on the cluster galac- 
tocentric distance for mass-segregated and unsegregated 
clusters. At a given i?g, the amount of impulsive mass 
loss is the same for both segregated and unsegregated 
clusters but, as anticipated by the analytical calculations 
presented earlier, the destructive effect of this mass loss 
is strongly augmented by mass segregation. 

Fig. m shows that all of the unsegregated clusters sur- 
vive the early mass loss due to stellar evolution and even- 
tually dissolve as a result of the evaporation of stars due 
to two-body relaxation. On the other hand, early impul- 
sive mass loss leads to the dissolution within just a few 
dynamical times of all of the initially mass-segregated 
clusters. Only the dissolution time of the cluster closest 
to the Galactic center (for which the dynamical time is 
short and the impulsive mass loss is negligible) is similar 
to that of the corresponding unsegregated cluster. 

As already noted in a number of previous studies (see 
e.g. Vesperini & Heggie 1997, Portegies Zwart et al. 
1998, Baumgardt & Makino 2003), the dissolution time 
of tidally truncated unsegregated clusters increases lin- 
early with Rg. This dependence is evident from the lower 
panel of Fig|5l in the upper panel it is incorporated in the 
dynamical time, tud, used to normalize the Tdiss (as we 
pointed out above in this section, tud oc Rg)- The disso- 
lution of the unsegregated clusters is driven principally 
by two-body relaxation, and for clusters more massive 
than those we have studied in our simulations, the ratio 
Tdiss /itid increases with the initial number of stars A^ as 
~ {N/logNY with X - 0.75 - 0.8 (see Baumgardt & 
Makino 2003). 

For initially mass-segregated systems the dependence 
of Tdiss on Rg is determined by the initial degree of mass 
segregation, the ratio R/Rt and the amount of impulsive 
mass loss. As noted above, for example, for clusters like 
S-Rgl with small dynamical times, the amount of impul- 
sive mass loss is negligible and the dissolution time of a 
mass-segregated cluster is similar to that of an unsegre- 
gated cluster. On the other hand, for very large amounts 
of impulsive mass loss — larger than those in the systems 
we have studied — the dissolution time must converge to 
a value of ^ tud- It is therefore not possible to derive a 
single general power-law to fit the scaling of Tdiss with 

Rg. 
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Fig. 4. — {Upper panel) Dissolution time, T^ias, (in units of 
cluster dynamical time, tfi^) versus galactocentric distance, Rg, 
for the segregated (red filled dots-5 runs) and unsegregated clus- 
ters (black open dots-f/ runs). The green open triangle and 
the blue asterisk-dot show, respectively, the dissolution time for 
the moderate-segregation {M-Rgl8) and low-segregation (L-Rgl8) 
runs. The black filled diamond-dot and the open orange diamond- 
dot show, respectively, the dissolution times for S-Rgl8R075 and 
U-Rgl8R 075; the purple cross and the purple square-dot (al- 
most overlapping) show, respectively, the dissolution times for S- 
Rgl8R05 and U-Rgl8R05 (for a consistent comparison of all the 
models at Rg = 18 kpc, the dynamical time used to normalize the 
dissolution times of S-R5iSfl075, U-Rgl8R 075, S-Rgl8R05 and U- 
Rgl8R05 is the same used for the other models at Rg = 18 kpc 
and equal to the dynamical time of the tidally truncated systems). 
The small red filled dot, the small green filled triangle (overlapping 
with the large triangle) and the small blue asterisk-dot show the 
dissolution times for S-Rgl8, M-Rgl8 and L-Rgl8 repeated with a 
larger number of particles {N = 70K, = 103K and N = 120K re- 
spectively). {Lower panel) Dissolution time (in Myr) versus galac- 
tocentric distance (symbols as in the upper panel). 



For the specific set of standard segregated S runs, and 
limiting the fit to systems with Rg = 4, 18, 40, the 
dissolution time is much shorter than that of the unseg- 
regated systems and has a weaker dependence on Rg, 

^diss ^ Rg 

The dissolution times of the mass-segregated clusters 
undergoing rapid dissolution do not depend on N , as 



Fig. 5. — Time evolution of the 10, 25, 50, 75, 90 percent La- 
grangian radii for S-Rgl8 (red solid line) and U-Rgl8 (black dashed 
line). 



shown by the results of simulations starting with the 
same degree of segregation but larger values of N (see 
Fig.il). 




2000 4000 6000 8000 10000 12000 
time {Myr} 

Fig. 6.— Time evolution of the 10, 25, 50, 75, 90 percent 
Lagrangian radii for S-Rgl8R05 (red solid line) and U-Rgl8R05 
(black dashed line) at Rg = 18 kpc. 



The difference in response of segregated and unseg- 
regated clusters to the same amount of impulsive mass 
loss is further illustrated by Fig. O which shows the time 
evolution of the 10, 25, 50, 75, 90 percent Lagrangian 
radii for S-Rgl8 and U-Rgl8. Both clusters initially lose 
the same fraction of their initial mass and expand in 
response to this mass loss; however the preferential re- 
moval of mass from the innermost regions in the segre- 
gated cluster, leads to a much stronger expansion and to 
rapid cluster dissolution. 

The longer dissolution times L-Rgl8 and M-Rgl8 plot- 
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ted in Fig. 2] show that, as the initial degree of mass 
segregation decreases, so does the disruptive effect of the 
early mass loss due to SN ejecta. 

The long-lived globular clusters we observe today 
might have had very low degrees of initial mass segre- 
gation which allowed them to survive. Alternatively, a 
cluster with a higher degree of initial mass segregation 
can survive if its initial size is significantly smaller than 
its tidal radius. In this case, most cluster stars remain 
within the tidal boundary during the early expansion, 
and the number of stars lost during this phase is smaller. 

To illustrate this point, we have repeated the S-Rgl8 
and U-Rgl8 runs, but setting the ratio of the initial 
cluster size, R, to the tidal radius Rt equal to 0.75 (5"- 
Rgl8R075 and U-Rgl8R075) and 0.5 {S-Rgl8R05 and 
U-Rgl8R05) . The dissolution times for these systems 
are plotted in Fig. [4] and show that the lifetime of a 
mass-segregated cluster with the same degree of initial 
segregation increases as the ratio R/Rt decreases. 

In particular, the lifetimes of S-Rgl8R05 and U- 
Rgl8R05 are very similar, as the early mass loss disrup- 
tive effects for the initially mass-segregated system are 
strongly suppressed by the fact that the cluster is ex- 
panding well within its tidal radius and not immediately 
losing stars as it expands. Since the system is initially 
denser and has a shorter dynamical time, the amount of 
impulsive mass loss is smaller than for the tidally trun- 
cated system; this effect also contributes, although to a 
lesser extent, to the extension of the cluster lifetime. 

The evolution of the 10, 25, 50, 75, 90 percent La- 
grangian radii for S-Rgl8R05 and U-Rgl8R05 are shown 
in Fig. [51 Although the lifetimes of these two systems 
are similar. Fig. [6] clearly shows important differences in 
the structural evolution of the two clusters: the mass- 
segregated cluster undergoes a stronger initial expansion 
which, although it does not lead to its prompt disso- 
lution, significantly delays the cluster core collapse and 
keeps the cluster concentration lower than that of the un- 
segregated cluster. This is further illustrated by Fig. [7] 
which shows the time evolution of the cluster concentra- 
tion as measured by the ratio of the core and half-mass 
radii, rc/ru, for S-Rgl8R05 and U-Rgl8R05. 

3.3. Results: evolution of the stellar mass function 

An important issue related to a cluster mass loss is 
the evolution of the stellar MF. It has been shown in a 
number of investigations of initially unsegregated clusters 
(see e.g. Vesperini & Heggie 1997, Baumgardt & Makino 
2003) that mass loss due to two-body relaxation leads 
to the preferential escape of low-mass stars and to the 
flattening of the stellar MF. 

In particular, for long-lived clusters losing mass due 
to two-body relaxation the extent of the flattening of 
the slope of the stellar MF is strictly correlated to the 
fraction of the initial cluster mass left (see e.g. Fig. 16 
in Vesperini & Heggie 1997 and Eq. 14 in Baumgardt & 
Makino 2003). 

For long-lived initially unsegregated clusters, the dy- 
namical mass loss timescale is closely related to the two- 
body relaxation timescale and is, in general, much longer 
than the stellar evolution timescale of the most massive 
stars. In this case, the range of stellar masses involved in 
segregation and evaporation processes after a few billion 
years includes the most massive dark remnants (1.2-1.4 
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2000 4000 6000 8000 10000 12000 
time (Myr) 

Fig. 7. — Time evolution of the ratio of the core to the half-mass 
radius, r^/r^, for S-Rgl8R05 (red solid line), U-Rgl8R05 (black 
dashed line). 



m,(7)), the more numerous white dwarfs, and the remain- 
ing main sequence stars (m < Im©). The most mas- 
sive main sequence stars usually included in the calcu- 
lation of the slope of the mass function fall in the high- 
mass end of this range, and hence tend to sink to the 
cluster center, while low-mass main sequence stars (e.g. 
m 0.1 — 0.2 TO©) preferentially escape from the cluster, 
flattening the MF. 

For an initially mass-segregated cluster, on the other 
hand, the range of stellar masses affected by segre- 
gation is much broader. The most massive stars in 
this broad range preferentially populate the innermost 
regions, while the long-lived main-sequence stars with 
(m < Im©) all have very similar spatial distributions. 

This point is illustrated in Fig. [H which shows the 
initial mean distance of stars from the cluster center as 
a function of the mean stellar mass in a number of mass 
bins for the S and U runs. The same information later in 
the evolution of the U-Rgl8R05, when a significant de- 
gree of mass segregation has been reached and the clus- 
ter hosts stars with a narrower range of masses, is also 
shown. 

Fig. [8] shows that the initial degree of segregation for 
stars with masses in the range 0.1 < Tn/m© < 0.5 (here- 
after we focus on the evolution of the MF of main se- 
quence stars with initial masses in this range) is negligi- 
ble and similar in clusters with and without primordial 
segregation. As discussed above, primordial mass seg- 
regation affects mainly the most massive stars initially 
present in the cluster. The segregation profile attained 
later in the evolution by an unsegregated cluster shown 
in Fig. [5] clearly illustrates the differences with the initial 
segregation profile of a cluster with primordial segrega- 
tion. This is an important point in understanding the 
implications of primordial mass segregation for the evo- 
lution of the MF. 

As the results presented in the previous sections show, 
the fate of initially mass-segregated clusters and the pro- 
cesses driving their dissolution depend on the degree of 
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U runs, t= Gyr 
S runs, t= Gyr 
S-Rg18R05, t^lGyr 
U-Rg18R05, t^lGyr 
U-Rg18R05, t-10 Gyr 
BDK08, t-0 Gyr 
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1.0 2.0 
<m> 
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10.0 20.0 
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time (Myr) 



10000 12000 



Fig. 8. — Mean radial distance from the center of the cluster of 
stars (normalized to the mean radial distance of stars in the first 
bin) binned in mass versus the mean mass of stars in each bin. The 
different lines show (from the top to the bottom line) the U runs 
at i = (dashed line-open dots), the S runs (solid line-filled dots) 
at t = 0, S-Rgl8R05 at t = f Gyr (red solid line-asterisk dots) 
U-Rgl8R05 at t = f Gyr (black dashed line-open square dots) and 
U-Rgl8R05 (black dashed line-crosses) at t = fO Gyr (the lines for 
the models at 4 = f Gyr are almost completely overlapped). The 
blue dotted line- / open triangles shows the initial segregation of one 
of the models with initial segregation (model No. 3 in their Table 
2) studied by Baumgardt, De Marchi & Kroupa (2008) (BDK08). 



initial mass segregation and on a cluster's initial size rela- 
tive to its tidal radius. Hereafter we focus on the S-Rgl8, 
S-Rgl8R05, U-Rgl8, C/-%i,Si?05 runs, to explore the dif- 
ference between the MF evolution in systems with and 
without primordial mass segregation. 

We fit the mass function of main sequence stars with 
masses in the range (0.1 < m/wQ < 0.5) with a power 
law dA^(TO) = Am~°'dm\ the upper panel of Fig. [9] 
shows the time evolution of a for S-Rgl8 and U-Rgl8. 
As discussed in the previous section, the initially mass- 
segregated system rapidly dissolves as its outer layers ex- 
pand in response to mass loss due to SN ejecta. As the 
system quickly loses mass, the MF slope also rapidly flat- 
tens; the evolution, mass loss, and MF flattening occur 
on a much longer timescale for the unsegregated system. 

The lower panel in Fig[5]shows the MF slope versus the 
fraction of the initial mass left in the cluster, clearly illus- 
trating an important difference between mass-segregated 
and unsegregated clusters: for a given amount of mass 
loss, the degree of MF flattening is significantly larger 
for the initially unsegregated cluster. This is due to the 
differences in the range spanned by stellar masses during 
the cluster evolution. For the slowly evolving initially 
unsegregated cluster, the narrower range of star masses 
present in the cluster during much of its evolution results 
in greater preferential loss of stars at the low-mass end 
of the range 0.1 < m/mQ < 0.5, while stars at the high- 
mass end of this range sink to the cluster center. This 
is not the case for the initially mass-segregated system, 
in which rapid dissolution ensures that stars much more 
massive than 0.5toq are present during most of the evo- 
lution, preventing the development of strong differences 
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Fig. 9. — {Upper panel ) Time evolution of the slope of the 
mass function, a, for main sequence stars with 0.1 < m/niQ < 
0.5 for S-Rgl8 (red solid line), U-Rgl8 (black dot-dashed Une), 
S-Rgl8R05 (red dashed Une) and U-Rgl8R05 (black dotted Une). 
{Lower panel) a versus the fraction of the cluster initial mass left 
in a cluster (different lines refer to the same runs as in the upper 
panel) . 



in the mass loss rate and the spatial segregation of stars 
with masses in the range 0.1 < m/mQ < 0.5. 

The time evolution of a and the dependence of a on 
the fraction of mass remaining in the cluster for the long- 
lived S-Rgl8R05 and U-Rgl8R05 systems are also shown 
in Fig.O The evolution of a is very similar for these two 
systems. The initial segregation profiles for U-Rgl8R05 
and S-Rgl8R05 in Fig. [8] clearly show that the degree 
of mass segregation for stars with masses in the range 
0.1 < m/mQ < 0.5 is very similar in these two systems. 
If a mass-segregated cluster survives its early expansion, 
as in the case of S-Rgl8R05, it will enter the evolution- 
ary phase driven by two-body relaxation (after the most 
massive stars have evolved off the main sequence) with 
a segregation profile similar to that of an initially un- 
segregated cluster (see the mass segregation profiles of 
S-Rgl8R05 and U-Rgl8R05 at t = 1 Gyr in FigE]), and 
the MFs of both systems will evolve similarly. 
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For clusters initially underfilling their Roche lobes, our 
simulations show that a given amount of mass loss leads 
to a stronger MF flattening than for tidally truncated 
unsegregated clusters. This is a consequence of the slower 
mass loss rate of clusters initially lying inside their tidal 
radii; it takes longer for these clusters to lose a given 
amount of mass, and during this additional time a higher 
degree of mass segregation for stars with masses 0.1 < 
m/m© < 0.5 is reached. 

4. DISCUSSION AND CONCLUSIONS 

The results of our study show that initial mass segre- 
gation can significantly affect the dynamical evolution of 
star clusters. 

For clusters initially filling their Roche lobes, our sim- 
ulations show that mass-segregated clusters can quickly 
dissolve as a result of the early expansion triggered by 
mass loss due to stellar evolution, unless the initial degree 
of mass segregation is low. We have studied the evolution 
of mass-segregated and unsegregated clusters undergoing 
the same amount of early impulsive mass loss and find 
that, for segregated clusters, the preferential removal of 
this mass from the innermost regions strongly augments 
the strength of the early expansion and leads to rapid 
cluster dissolution. 

We have also explored the evolution of clusters initially 
underfilling their Roche lobes, showing that in this case 
mass-segregated clusters can survive the early expansion 
phase; a mass-segregated cluster still undergoes a strong 
expansion in this case, but much of this expansion occurs 
within its tidal radius without significant loss of stars. 
The subsequent cluster evolution, lifetime, and mass loss 
rate do not differ significantly from those of an initially 
unsegregated cluster. 

The stronger expansion of the segregated cluster, while 
not leading to rapid dissolution, does drive the cluster to- 
ward a less concentrated structure, significantly delaying 
core collapse. The less concentrated structure of an ini- 
tially segregated cluster throughout most of the evolution 
is the main difference between long-lived segregated and 
unsegregated clusters. If clusters still surviving today 
were initially segregated, this is one of the possible evo- 
lutionary routes they might have followed. Another pos- 
sible scenario has been studied by D'Ercole et al. (2008), 
who demonstrate the role of primordial segregation in the 
evolution of clusters which form a second generation of 
stars in the innermost regions of a first-generation popu- 
lation. The early expansion in this case plays a key role in 
causing the loss of a large fraction of the first-generation 
stars, while most of the centrally concentrated second 
generation remains in the cluster, resulting in a system 
containing a mix of first and second generation stars that 
eventually survives until the present day. 

We have also studied the differences in the evolution 
of the stellar MF for segregated and unsegregated clus- 
ters, focusing on the evolution of the MF of low-mass 
(0.1 < m/rriQ < 0.5) main sequence stars. Our simula- 
tions show that the MF of a rapidly dissolving segregated 
cluster flattens quickly as the system loses stars. How- 
ever the extent of the MF flattening induced by a given 
amount of mass loss for a tidally truncated segregated 
cluster (our S runs) is smaller than for an initially un- 
segregated system (our U runs). 

The difference is a consequence of the rapid dissolution 



of segregated clusters. We have shown that the initial 
degree of segregation adopted in our simulations mostly 
affects massive stars with m > 5 — IOtoq, while low- 
mass stars with m < Im©, initially have a similar spa- 
tial distribution to stars in unsegregated clusters. Seg- 
regated clusters initially filling their Roche lobe quickly 
lose mass before two-body relaxation can produce a sig- 
nificant difference in the spatial distribution (and in the 
mass loss rate) of stars with masses m < Im©. Unsegre- 
gated clusters, on the other hand, lose stars more slowly 
and a large fraction of their initial mass is lost late in 
the evolution, when the most massive stars in the range 
0.1 < m/rriQ < 0.5 are significantly more concentrated 
than stars of lower masses. The differences between the 
iS" and U systems in the degree of MF flattening for a 
given amount of mass loss result from differences in the 
processes driving the cluster dissolution and, more im- 
portantly, from differences in the dissolution timescales 
associated with these processes. 

Segregated clusters initially well conflned within their 
Roche lobe may survive their early expansion and, in 
this case, we have shown that both segregated and un- 
segregated clusters display very similar MF evolution. As 
pointed out above, the main fingerprint of initial segre- 
gation might be looser cluster structure during a larger 
fraction of the cluster lifetime. However, we emphasize 
that, unless some other mechanism continues to drive 
the cluster expansion, long-lived segregated clusters also 
eventually undergo core collapse; therefore it can not be 
ruled out that high-concentration and post-core collapse 
clusters might have been initially segregated ^. 

It has recently been suggested by Baumgardt et al. 
(2008) that the low-concentration clusters with strongly 
fiattened MFs found in the observational study by De 
Marchi et al. (2007) can result only from the dynamical 
evolution of initially segregated clusters (see also Marks 
et al. 2008 for the possible role of primordial gas expul- 
sion). In particular, the simulations of Baumgardt et al. 
show that a cluster must initially be significantly mass 
segregated and fill its Roche lobe in order to evolve into 
a system with properties resembling those of the loose, 
fiat-MF clusters observed by De Marchi et al. The ini- 
tial degree of segregation adopted by Baumgardt et al. 
is significantly larger than that of our simulations, as 
shown in FigO However, the Baumgardt et al. study 
focuses on the effects of two-body relaxation, and does 
not include the early mass loss due to stellar evolution 
which, as shown by our simulations, can play a crucial 
role in cluster evolution. As discussed above, initially 
tidally limited segregated clusters may dissolve rapidly 
due to the expansion triggered by early mass loss, even 
when starting with a degree of initial segregation much 
smaller than that adopted by Baumgardt et al. 

In a future investigation we will focus on detailed com- 
parisons between our simulations and observational data, 
with the goal of assessing which, if any, observed cluster 
properties might represent unique signatures of initial 
mass segregation. 

1 See e.g. Mackey et al. (2007, 2008) for a study of the early 
expansion of segregated clusters driven by stellar evolution mass 
loss and the subsequent heating from a population of stellar black 
holes, and the role these processes might play in determining the 
radius-age trend observed for massive clusters in the Magellanic 
Clouds. 
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